Topological Winding and Unwinding in Metastable Bose-Einstein Condensates 
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Topological winding and unwinding in a quasi-one-dimensional metastable Bose-Einstein con- 
densate are shown to be manipulated by changing the strength of interaction or the frequency of 
rotation. Exact diagonalization analysis reveals that quasidegenerate states emerge spontaneously 
near the transition point, allowing a smooth crossover between topologically distinct states. On a 
mean-field level, the transition is accompanied by formation of grey solitons, or density notches, 
which serve as an experimental signature of this phenomenon. 
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Metastability of a physical system leads to a rich va- 
riety of quantum phases and transport properties that 
are not present in the ground state phase. An illus- 
trative example is superflow and phase slip in a nar- 
row superconducting channel [l| . Other examples include 
Feshbach molecules formed in high rotational states [2[ 
and metastable quantum phases in higher Bloch bands 
in an optical lattice [3(]. Recent experimental advances 
in cold atoms/molecules have made it possible to real- 
ize excited, metastable states which persist for a long 
time. These states provide an excellent medium in which 
to investigate fundamental aspects of condensed matter 
systems such as topological excitations and superfluid- 
ity @,S,i,0. 

It is widely believed that the angular momentum per 
particle in a weakly repulsive one-dimensional (ID) su- 
perfluid ring system @, 0] is quantized at T = and 
that there are discontinuous jumps between states hav- 
ing different values of the phase winding number. In this 
Letter, we point out that this applies only to the ground 
state; continuous transitions do in fact occur between 
metastable states of repulsive condensates. The underly- 
ing physics behind this phenomenon is the emergence of 
a dark or grey soliton train [8] which bifurcates from the 
plane-wave solution and carries a fraction of the quan- 
tized value of the angular momentum. 

Starting with mean-field theory for scalar bosons sub- 
ject to rotation, we proceed through progressively deeper 
levels of insight into the quantum many-body nature of 
this problem, making a link between semiclassical and 
quantum solitons in metastable states. We find that 
the phase slip, which allows a smooth crossover be- 
tween topologically distinct states, is caused by a quan- 
tum soliton. The latter consists of a linear superposi- 
tion of the rotationalfy-invariant many-body eigenstates 
of the Hamiltonian [9[. In both Boguliubov theory and 
quantum many-body theory the broken-symmetry soli- 
ton state is shown to be stable against perturbation. 

This phenomenon can be realized by hot atoms con- 



fined in fast-rotating circular waveguides or toroidal 
traps [10[. First, to obtain a metastable uniform con- 
densate one quickly stops the rotation and then lowers 
the temperature. Second, one adiabatically changes the 
angular frequency of the trap in the presence of a small 
arbitrary perturbation in the trapping potential. This 
causes atoms to adiabatically take the higher-energy path 
of a metastable soliton state, as we will show. Third, one 
stops the adiabatic change in the frequency at the cor- 
rect point to arrive at a different winding number. All of 
these processes can occur continuously. 

We consider a system of N bosonic atoms in a quasi- 
1D torus with radius R, under an external rotating drive 
with angular frequency 2^. The length, angular momen- 
tum, and energy are measured in units of R, h, and 
H 2 /{2mR 2 ) 1 respectively. The Hamiltonian is given by 
the Lieb-Liniger Hamiltonian in a rotating frame of ref- 
erence [111 Il2j. 

H = dO[^(-id e - fi) 2 ^ + gm^H 2 /% (1) 

where #id characterizes the strength of the s-wave inter- 
atomic collisions in ID [l3[ rescaled by H 2 /(2mR), 6 is 
the azimuthal angle, and the bosonic field operator sat- 
isfies periodic boundary conditions: ip(6) = ip(6 + 2tt). 
Since the Hamiltonian is periodic with respect to £7, the 
properties of the system are periodic in Q with period 
1 [14| , in direct analogy to the reduced Brillouin zone in 
a Bloch band @, EH . Without loss of generality we will 
henceforth restrict ourselves to Q G [0, 1). The Hamilto- 
nian is integrable via the boson-fermion mapping in the 
Tonks-Girardeau (TG) limit ^ D » iV 0, and via the 
Bethe ansatz in the weak-inter ation limit gir>N < O(l) 
as well as intermediate- interaction regimes. 

We first show how continuous changes in the angular 
momentum occur in the weak-interaction regime for so- 
lutions of the Gross-Pitaevskii equation (GPE) [(—ido — 
ty 2 +gmN\ilj(0)\ 2 }tJj(0) = ^(0), where ^ is the order pa- 
rameter normalized to unity and ip = Arg(^) is its phase. 
The single- valuedness of the wave function requires ip(0 + 
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FIG. 1: Amplitude (solid curves with the left reference), and phase (dotted curves with the right reference) of metastable 
states of the GPE for gmN = 0.6tv. Uniform solutions with different values of the phase winding (i) J = 1, and (vi) J = are 
smoothly connected through the broken- symmetry grey soliton (ii)-(v) with a self-induced phase slip at ft = 0.5. 



2tt) = cp(0) + 2?r J, where J G {0, ±1, ±2, . . .} is the topo- 
logical winding number. Stationary solutions of the GPE 
for #id > are either plane- wave states i/)(0) 



r)E)]; the phase 



sn /i sn /2/ sn ; there are j density 



or a grey soliton train [8[ whose amplitude and phase 
are given by \^{6)\ = A[l + 77 dn 2 (jK(6 - o )/tt, /c)] 1/2 , 
and <p(0) =^6 + B U (&jK(0 - 6> )/tt, k), respectively. 
Here the amplitude A = ^K/[2tt(K 
pre-factor B = (S/jK)y/g { 
notches in the soliton train; 77 = —2j 2 K 2 /g sn G [—1,0] 
characterizes the depth of each density notch; k G [0, 1] is 
the elliptic modulus; K(k), E(k), fc), are elliptic 

integrals of the first, second, and third kinds; and dn(?i, k) 
is the Jacobi dn function. The degeneracy parameter 0q 
indicates that the soliton solutions are broken-symmetry 
states. We define / sn = irg 1D N/2 - 2j 2 K 2 + 2j 2 KE, 
gsn = /sn + 2j 2 K 2 , h sn = f sn + 2k 2 j 2 K 2 , and S = 1 
for < ft < 0.5, 5 = -1 for 0.5 < ft < 1. Then 
£ = —2(kjK) 2 /f sri < , and only when soliton solutions 
exist [17] is fc 7^ 0. In the limit rj —> —1, f sn approaches 
zero, and g sn and h sn both approach finite positive values; 
consequently, the wave function approaches the Jacobi sn 
function, which corresponds to a dark soliton with a tt- 
phase jump at 0q. In the opposite limit r/ — > 0, both the 
amplitude and phase approach those of the plane-wave 
solution with the same phase winding J. These limiting 
behaviours indicate that continuous change in topology 
of the condensate wave function is possible, as illustrated 
in Fig. [1] (i)-(vi). Henceforth, we consider the single soli- 
ton j = 1 for simplicity, but our discussion holds for 
arbitrary soliton trains j > 1. 

Bifurcation of the soliton train from the plane wave 
constitutes a second-order quantum phase transition with 
respect to #id and/or £1. Figure [21(a) shows the energy 
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(«- J) 2 +£i D A7(4^), 

<?i D A7(27t) + [3KE - K 2 (2 - k 2 )] /tt 2 



+4K 2 [3E 2 -2(2 - k 2 )KE+K 2 (l-k 2 )]/(37r 3 g 1B N). 

This kind of bifurcation does not occur from the ground- 
state energy. However, for metastable states a bifurca- 
tion can occur between the plane- wave state and the soli- 



ton state with the same winding number J. After bi- 
furcation, the soliton energy Ej* becomes larger than 



^ pw) . Furthermore, at ft = 0.5, ^ aui; and E{* UL) are de- 
generate with a ±7r-phase jump in the condensate wave 
function. The derivatives of the energies dE { J ol) /dn and 
QE& w )/dQ have a kink at the boundary as can be ver- 
ified analytically. This identifies the second-order quan- 
tum phase transition [18], which occurs along a curve in 
the ft-#iD plane. 

Figures [U (i)-(vi) illustrate a continuous change in the 
topology along a higher-energy, soliton path shown in 
Fig. [2] (a) with white arrows. Following this path in 
Fig. [TJ as ft increases starting from (i) the plane wave 
with J = 1, (ii) solitons start to form past a critical point 
ft cr . (iii) The density notch deepens for ft cr < ft < 0.5. 
At ft = 0.5 it forms a node, the phase of the soliton 
jumps by 7r, and the energies of the solitons with phase 
winding number 1 and are degenerate, (iv), (v) The 
soliton with phase winding J = deforms continuously 
as ft increases, (vi) Finally, the state goes back to the 
plane- wave state with phase winding J = [l9l ]. 

The angular momentum L/N = J dOtp* (—ide)i(j of the 
metastable states changes continuously along the soli- 
ton path. For the plane wave state, Lj PW ^/7V = J 
is quantized; in contrast, for the soliton Lj 01 ^ /N = 
ft + S\/2f sn g sn h sn / (gir>N7r 2 ) is non-integer, as shown 
in Fig. [2Kb). Thus a continuous change of angular mo- 
mentum is possible for ID Bose systems by taking the 
metastable states with energy slightly higher than that 
of the ground state. 

We next investigate the stability of the metastable 
states using Boguliubov theory 20|, 21], and identify the 
curve in the ft-giD plane where the soliton solutions bifur- 
cate from the plane- wave solutions. A stationary solution 
ip of the GPE subject to a small perturbation S evolves 
in time as $(t) = e~^[^ + E n (^e"' Ari£ + Sv*e iX *% 
where (u n , v n ) and A n are eigenstates and eigenvalues of 
the Boguliubov-de Gennes equations (BdGE), and n de- 
notes the index of the eigenvalues. 

For the plane-wave state with phase winding J, the 
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FIG. 2: (a) Energy difference between the metastable plane- 
wave and soliton states. The soliton solutions exist in the area 
sorrounded by two phase boundaries (white dotted curves), 
(b) Corresponding angular momentum. The soliton solutions 
make it possible to smoothly connect quantized integer values 
of average angular momentum. 



negative, monotonically increases for ft G [0,ft cr ), and 
crosses zero at (giD^V)cr/(27r) - 2(ft cr - J) 2 + 1/2 = 0. 
Thus the metastable state ip = e ze /^/27r is thermody- 
namically unstable. The plane-wave limit of the soliton 
solutions r] — > occurs when (gmN, approach their 
critical values from above. We also obtained the eigen- 
values of the BdGE when ip is taken as a soliton. The 
eigenstates obtained from a soliton involve the Nambu- 



Goldstone mode A 



(J=l,sol) 
NG 



i.e., the zero-energy rotation 



mode, associated with the rotational symmetry breaking 
of the soliton solution. At the critical values of (ginN, ft), 



A 



(J=l,pw) 



Avn ±,bUi; and other eigenvalues A^ J ' so1 ^ are 



(J=l,sol) 
1 — /X NG 

nearly degenerate with A^ J ' pw ^. There is no negative 
Boguliubov mode in the soliton regime. Thus the soli- 
ton state is linearly stable. 

Finally, we investigate how continuous change in the 
angular momentum per particle observed in the mean 
field theory is described in terms of quantum many-body 
theory. Lieb derived two kinds of excitation branches in 
the thermodynamic limit, "Type I" and "Type II" exci- 
tation branches of Eq. (pp) with ft = [11]. The Type II 
excitation was shown to be a dark soliton branch [22] . We 
also determined low-lying excitation energies in a finite- 
size waveguide in a wide range of groN via the Bethe 
ansatz [llHiH. 

In order to study the metastable states we need to 
investigate the 7V th order highly excited states. We di- 
agonalize Eq. (pp) in a basis which provides equivalent 
results to those obtained with the Bethe ansatz [26[ for 
the purposes of physical insight. The basis is taken as 
|n_i, no, ni, 77,2), subject to conditions rt\ = TV and 
Ini = L, where n\ is the number of atoms with single- 
particle angular momentum I and L G {—N, . . . , 2 N} is 
the total angular momentum. Figure [3](a) shows the en- 
ergies El = (L, N\H\L, N) of yrast states |£,iV), i.e., 
the lowest-energy state for a fixed value of L and TV. 
The index N will be dropped, as it is fixed. The curva- 
ture of the surface is independent of N but the density 



of states with respect to L/N increases as TV becomes 
larger. Three kinds of states appear in this energy land- 
scape: 

(i) Ground state: The ground state of the Hamiltonian 
is \L = 0) with energy E L=0 /N ~ ft 2 + gmN/(4,n) for 
ft G [0,0.5), and \L = N) with energy E L=N /N ~ (ft - 
l) 2 + #idN/(4tt) for ft G [0.5,1), respectively. These 
are in agreement with the mean-field ground states ip = 
l/y/27r and ip = e %e /y/27r, respectively. 

(ii) Metastable plane-wave states: Similarly, the many- 
body counterparts of the metastable plane- wave solutions 
of the GP equation are \L = N) for < ft < 0.5 and 
\L = 0) for 0.5 < ft < 1, respectively. 

(iii) Seeds of broken-symmetry states: All many- 
body eigenstates are rotationally invariant because they 
respect the symmetry of the Hamiltonian, i.e., there 
are no broken-symmetry eigenstates. However, as 
shown in Fig. [3] (b), the eigenvalues {El=q, . . . , Ejy} 
(densely packed blue curves) cross each other in be- 
tween the ground and metastable plane- wave states, and 
{El=\, . . . , i£jv-i} become larger than the metastable 
plane- wave state energies within a certain range of ft. 
The regime of this quasi-degenerate level crossing can be 
identified with the soliton regime given by the mean-field 
theory, as we will show. Furthermore, the envelope of the 
highest eigenvalues in the quasidegenerate regime coin- 
cides with i£j°o 1 in the limit TV — > 00. In the absence 
of interaction (#id = 0), all the levels with L G {0,7V} 
are degenerate for ft = 0.5 only. The degeneracy inher- 
ent in a level crossing permits the existence of a broken 
symmetry state; in our case, this is the soliton solution 
corresponding to the mean-field theory prediction. 

Figure [3] (c) shows the two-body correlation func- 
tion g^te-e') = (^^o)^^o , )^(o , )^(o)}/[(^ (0)^(0)} 

(ijj^ (6 f )ip(6 f ))], where the expectation value is taken with 
respect to each yrast state \L). The function gW is 
independent of ft and has a single peak. When L/N 
is an integer, g^ = 1; when L/N is half- integer, 
g^ strongly deviates from 1. The one-body density 
g^ = (ijj^ (0)ip(0)} = 2tt/N is a constant for any eigen- 
state, since the system is rotationally invariant. 

In order to force quantum solitons to appear in g^\ 
we add a symmetry-breaking perturbation of the form 

V = £Eze{-i,o,i,2}(4fi^ + hx 0> £ < 9ib- The total 
angular momentum L is no longer a good quantum num- 
ber, and the eigenvalue problem is thus given in a gen- 
eral form, (H + V)\^ n ) = E n \^ n ) where n G {0, 1, . . . } 
is the energy eigenvalue index. Employing the basis |L), 
where L G {—TV, 27V}, the average angular momentum 
per particle is shown in Fig. [U(d) for each eigenstate for 
e = 5 x 10 -3 . Outside the soliton regime, angular mo- 
menta remain integer values unaffected by the small per- 
turbation. In the soliton regime, on the other hand, by 
a superposition of eigensates \L) in the absence of the 
perturbation some of the angular momenta converge to 




FIG. 3: (a) The lowest energy of the Hamiltonian for each an- 
gular momentum subspace for N = 40, gm = 2tt x 7.5 x 10 -3 . 
(b) Enlargement of (a) near the critical point. Densely packed 
blue curves are spectra within L £ {0, N}. Otherwise the 
spectra are sparse, (c) Two-body correlation function of 
eigenstates for each angular momentum subspace. (d) Ex- 
pectation value of the angular momentum of each eigenstate, 
in the presence of symmetry breaking potential V as a func- 
tion of Q. 

a noninteger value, which indeed agrees with £j° v We 
have also calculated the one-body correlation function 
for all eigenstates in the presence of the symmetry 
breaking perturbation, and confirm that has a single 
density notch in the soliton regime with the depth close 
to the mean-field grey soliton. 

In conclusion, we found a denumerably infinite set of 
paths to connect plane- wave states via soliton trains in a 
metastable system of scalar bosons on a ring. Associated 
with this transition, the energy of the solitons bifurcates, 
and a continuous change in the angular momentum be- 
comes possible in the mean-field theory. We made a link 
between these mean-field results and the full quantum 
theory by showing that quasidegenerate energy levels are 
related to the formation of quantum solitons. 
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